Semiclassical theory for quantum spin models with ring exchange on the triangular 

lattice 
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A semiclassical theory of a quantum spin— 5 model with competing ring and Heisenberg exchange 
terms on the triangular lattice is obtained. A mechanism for the generation of Z2 vortices is 
exhibited. The vortices are shown to carry a nontrivial geometric phase for the order parameter 
when 2S is odd, leading to a difference between the quantum disordered ground states and low 
energy spectra for half odd integer and half even integer spin systems, and a topological degeneracy 
on surfaces with nontrivial cycles. A connection to dimer models is discussed. 



Since the original proposal 0, that frustrated quan- 
tum spin half systems on the triangular lattice may fail 
to exhibit any form of local order at zero temperature, 
many efforts have been made to characterize such spin 
liquid states Quantum spin (and bosonic) systems 
with global XY symmetry have been argued to possess 
disordered ground states with a proliferation of Z2 vor- 
tices Q, as have dimer models H B\ and Sp{N) mod- 
els for large— A'' Q- Exact diagonalizations §3 9] of the 
SU{2) invariant spin half model with ring exchange inter- 
actions on the triangular lattice indicate the lack of any 
local order or lattice symmetry breaking in the thermo- 
dynamic limit but an understanding of the microscopic 
origin of the disorder in such spin models is incomplete. 
In this paper a semiclassical approach to the low energy 
physics of the ring exchange model and similar SU (2) in- 
variant systems with noncoUinear order is presented 0, 
and the microscopic origin of Z2 vortices which disorder 
the system is elucidated. A geometric phase associated 
with the vortices, arising from a topological term in the 
effective action for spin 5, is argued to give rise to basic 
distinctions between the low energy spectra for the cases 
when 2S is odd or even. 

The Hamiltonian of the ring exchange model for spin 
half is 

H = J2 ^ (S, • S,) + J4 ^[4(S, • Sj)(Sfe • SO (1) 

<ij> ijkl 

+4(S,-S,)(S, •Sfe)-4(S,.Sfc)(S, -S,)] 

where < ij > denotes that the sites are nearest neigh- 
bours and in the second term the sum is over all elemen- 
tary plaquettes bounded by sites ijkl on the triangular 
lattice. These interactions have been used to describe 
He^ monolayers on graphite ^3 ^^'^ generated in 
order f^/U^ in the t/U expansion of the half filled Hub- 
bard model. Since each spin is shared by twelve rhombi 
the spins are highly frustrated. This model will be stud- 
ied for large S; to keep the ratio J4/ J2 fixed in this limit 
J4 must be scaled by S^. 

A low energy theory for this model is obtained by 
treating the effects of short range fluctuations at tree 



level (to order in 1/S), while maintaining the nonlinear- 
ity of the slowly varying fluctuations of the local clas- 
sical order and accounting for the Berry phases which 
arise in the path integral for the time evolution operator 
[ill . It was shown in that the classical limit of 
the ring exchange Hamiltonian (QJ has a four sublattice 
"tetrahedral" ground state on the triangular lattice for 
\J4/J2\ < 2. The spins on the vertices of an elementary 
rhombus A1A2A3A4 are given in terms of an orthonormal 
triad 6^,65, ec by 



{eA + es + ec), 
(e^ - bb - ec), 



_S_ 

— (-eA + es - ec), 
-^(-e^ - es -I- ec); 



(2) 



the spins point to the vertices of a regular tetrahedron. 
This state is invariant under global 0(3) x 0(3)spi„, un- 
der which the rotation matrix ea.a^ [ol ~ A, B, C; a ~ x, 
y, 2), transforms as e = L^eR, with L, R ^ 0(3); the 
additional 0(3) symmetry is special to the tetrahedral 
state; planar states in general possess 0(2) x 0{3)spin 
symmetry while an arbitrary nonplanar spiral state pos- 
sesses only the spin symmetry. 

The path integral for the time evolution operator, 
Z{T) = ex\)(—iHT)\n), may be obtained using spin 
coherent states [IJI by insertion of N resolutions of the 
identity in the matrix element for the time evolution op- 
erator in the coherent state basis, taking the limit of infi- 
nite A^ assuming that only continuous paths on the unit 
sphere contribute, and is given by 



Z= i:)r2e-s:i'''*'f^s(t)giS^(fj) 



(3) 



where the path integral is over closed paths in time, 
since it is obtained from diagonal matrix elements in the 
spin coherent state basis, ^lm,n = Sm,n/S is the nor- 
malized spin at site Xm,n — mui + nu2,with ui^2 being 
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elementary lattice vectors, Hs is the Hamiltonian (which 
may be arranged as series in S~^) and a;(f2) is the solid 
angle subtended by the closed path traced out on the 
unit sphere by the time evolution of ft. The solid angle 
uj{n) does not have an explicitly rotationally invariant 
local expression in fl; it is the action due to a monopole 
of strength 25 at the center of the unit sphere on which 
Cl lies. 

In the semiclassical limit the triad e^i, e^, is allowed 
to vary spatially; the spin is 



2S 



+ OiS-'), (4) 



where we define the orthonormal vectors 



-^((-l)™e^,,„,„ + (-l)"eB,„.,„ + (5) 

[—1) ec,7n,n), 



(-1)™+" 1 



-(-1) 



m+n 



2 ^C,m,n ) ; 



(-1) 



7n-\-n 



V2 



The short range fluctuation vector L is normal to ei; thus 

L = ((-l)'"Li + (-1)"L2 + (-l)™+"i3)e2,™.„ (6) 
+ {L^s + (-l)"+"i±,2)e3,™,„, 

where the vector L has been separated into compo- 
nents which vary rapidly in space with different frequen- 
cies. The action St — Lu{Sm,n) — § Hsdt is expanded 
to quadratic order in L and Ve, higher powers of L being 
lower order in S. The Berry phase term is 



Suj{Sm,nlS) = ^^(ei 



(7) 



dfLm.n ■ dtei,m,n X ei,m,n + 0[l/S). 



Thus, the spin path integral is 

^{Gl,7n,n) iSd) (8) 

where Sd = a Z]m,n / ^^I^™," ' c^tei,m,« x ei,m,„ - 
i/s(e,Ve,L)]. Since L appears quadratically in the ac- 
tion, it can be integrated out to obtain an effective action 
for the fields eA,B,c', terms in Sd which vary with mo- 
menta near tt times the reciprocal lattice vectors average 
to zero. In fact the 0(3) x 0{3)spin symmetry determines 
the sigma model action upto constants; on the lattice it 
is 

^ -Lp,e.(x„,^„)|2+ (9) 



Sd ^ jdt 



5)] 



where the (bare) couplings g"^ and a-i'e related to J2, 
J4 as above. These coupling will however flow with the 
length scale; in the following the effective couplings will 
be referred to. 

The evaluation of the Berry phase X)™ n '^(^i,™.") 
a special case of the result of .14, ISj and is done as fol- 
lows. Define the vectors f (x„j_„, u) and g(xm^„, v), which 
interpolate continuously on the unit sphere on the great 
circles between ei^m,n and ea^m.n, and between es^m.n 
and the slowly varying vector ec,m,n, by 

f (xm^„, u) = ei^m^„ sinu -I- e^^ra^n cosM (10) 
g(x„,„, u) = es^m^n sinw -I- ec,m,n cosw. 

The Berry phase at a given site, the solid angle 
Suj(ei^rn,n), Can be written as the solid angle swept out 
by the orthogonal vector e^^rn,n, Suj{ei^rn,n), plus the 
solid angle swept out by ei^m,n in a frame comoving with 



Si^{ei,m,n) = Suj{ez,m,n) + 

Si>dt /■^^^^^^^!il^.f(x„,„,^.)x 



(11) 

c)f(x„t^n,M) 

du 



Similarly, the solid angle swept out by the vector 63 „ „ 
can be expressed in terms of the solid angle swept out 
by the slowly varying vector ec,m,n orthogonal to e^^m.n 
using the interpolating vector g to obtain for the solid an- 
gle for ei_„i^„ an expression in terms of the slowly varying 
triad r c; 



Suj{ei,m,n) = Suj{ec^ra,n) + 



(12) 



S jdt j\ 
9f(x. 



c'g(x„i,„,u) ag(x„,„,it) 
•g(x,„,„,u)x 



dt 



du 



u) 9f(x„i_„,-u) 



ot ou 
and integrating over u, the Berry phase is 



Suj{ei^rn,n) — Suj{ec,m,n) — 
.rn,n 

dt 



(13) 



S <p dt-{eA, 



65, 



dsA 
di ' 



As shown in 'l4'| this quantity is a topological invariant; 
its variation under an arbitrarily small change in the triad 
e vanishes. If h.{eA,B.c} is defined by 



A{e^^s,c} = exp { - [w(ea, 



j> dt-{eA,m,n 



de 



B ,7n.n 



dt 



- Gb.,, 



deA,mM) 

dt 



(14) 

} 



then it can be seen that [lj| A{eA,B,c} G '^i{SO{3)), 
where tti (5*0(3)) is the fundamental group Z2 of 50(3). 
The spin path integral thus takes the form 

ZiT) = J DeA,B^cYl[HeA,BM^rn,n)}?^exp{iSdll5) 
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Consider in a given spin configuration entering in the 
path integral the spins in a line along a lattice vector. 
Labelling the time history of each spin along the line by 
^{^A,B,c} for that site, a change of A{eA^B,c} from one 
spin to the next requires a breakdown of continuity of 
the triad eA,B,c in spacetime. Analogously, introducing 
the spatial Z2 character along the chosen line at differ- 
ent time slices it can be seen that the breakdown of the 
triad occurs at a point in time. The point in spacetime 
where the breakdown of continuity occurs is the location 
of a (1 + 1 dimensional) instanton. When 25 is odd, the 
term Y[m,n[^{^A,B,G{'^m,n)}]'^^ appearing in the path in- 
tegral associates a phase tt to the rotation of the local 
triad about an axis in configuration space when 

the triad is transported around a loop enclosing the Z2 
instanton. In 2 -|- 1 dimensions, by continuing the instan- 
ton for the line in transverse directions, it is seen that 
Z2 vortices are present and that the phase holonomy of 
the local triad around any loop enclosing Uy vortices is 
2'KnvS. Since the triads must describe closed paths in 
configuration space in the coherent state path integral 
formulation, the vortex loops must be closed in the ab- 
sence of sources. 

Thus when 25 is even the effective theory is that of 
the 0(3) principal chiral model, which bears many simi- 
larities to the 0(3) nonlinear sigma model [IJj. For suf- 
ficiently small coupling go the principal chiral model de- 
scribes long wavelength fluctuations about the tetrahe- 
drally ordered state. In the ordered state the energy of a 
vortex loop is proportional to its length and to the spin 
stiffness. When go exceeds a critical value gc the spin 
stiffness vanishes, a gap opens up for local excitations 
and the vortex loops proliferate; their effects on the dis- 
ordered state depend also on the phase factors associated 
with them. 

When the system is placed on a two dimensional sur- 
face of genus g (e.g. by choosing appropriate bound- 
ary conditions) in the disordered phase the ground 
state degeneracy is determined by the Z2 character 
^m,n[^{^A,B,c{^m,n)}]^^ ■ Order parameter configura- 
tions which have a vortex which winds along a nontrivial 
cycle of the surface are not continuously deformable to 
those which do not contain the vortex; thus the classical 
configurations fall into 2g distinct sectors. Tunnelling 
between sectors can occur when a vortex pierces the sur- 
face. Consider such a Z2 vortex tunnelling between two 
such sectors through a plaquette; more complicated con- 
figurations of the vortex can be treated by decomposing 
the vortex into a sum of closed (contractible) loops and 
the above . Translation of the vortex by a lattice site 
in a transverse direction to the vortex leads to a phase 
change 27r5' in the path integral due to the phase for 
the triads associated with the vortex, but no change in 
energy. When 25 is odd there is a cancellation in the 
path integral between these configurations differing by 
the translation of the vortex by a lattice vector, and thus 



vortices are excluded from the system. This leads to 
a twofold degeneracy associated with each cycle and an 
overall degeneracy of the ground state of 2^^ when 25 is 
odd. 

If it is assumed that the disordered phase for g > gc 
can be continued to large g without a phase transition, 
the action can be studied in an expansion in g~^, as done 
in for the 1-1-1 dimensional 0(3) sigma model. In 
the large g limit the local kinetic energy of the triad at 
each site dominates the intersite couplings in the Hamil- 
tonian; i.e. this is a tight binding limit. It can be shown, 
for instance by computing in local coordinates on the 
configuration space of the triad 0| the Z2 character 
[A{eA,B,ci^m,n)}]'^^ at each site that it is the action due 
to a monopole of strength 2S{mod)2 at the origin of the 
solid sphere with opposite points identified (which is the 
configuration space of the local triad). In contrast to the 
case of a vector order parameter associated with antifer- 
romagnets the monopole charge is Z2 valued. As in the 
case of SU{2) spins the effect of the monopole is to de- 
termine the quantization of the symmetric top formed by 
the dynamics of the triad eA,B,c to be integer or half odd 
integer, according as 25 is even or odd. The eigenstates 
of a symmetric rotor are (2 J-|-l)^ fold degenerate; includ- 
ing the half odd integer sector they form the vector and 
spinor representations of 50(4). Thus when 25 is odd 
the eigenstates of the free rotor are classified by left and 
right angular momenta ( J^,, Jr) with both J^, Jfj being 
equal integers. The ground state is a nondegenerate sin- 
glet formed by the mixing with neighbouring sites by the 
hopping term in the tight binding model, and the lowest 
excited states carry quantum numbers {Jl, Jr) = (1, 1) 
and are ninefold degenerate. 

When 25 is odd the spectrum of the half odd integer 
rotor at each site has a fourfold degenerate ground state 
with [JliJr) = (515); restricting the Hilbert space at 
each site to these four states (a valid approximation for 
large g) and coupling adjacent sites through the hopping 
term leads to an effective Hamiltonian 

He = -\ ([^™™«^^m'n']'+ (16) 

<mn.7n' n' > 

-TrM'^ M ' ,) 

mn m n } 

where M is a generator of SO {A). Defining spin ^ oper- 
ators Sa, Ss by = \{S\ + S'b), M^^ = l^^'^iiS^ " 
5^) He = -\E<^J>i^A^ ' Saj){Sb. ' Sfi,). The four 
spin interaction term in He favours the formation of a 
singlet between neighbours, and the nearest neighbour 
dimer limit of this model reduces to that studied in |^ Q 
and is disordered (also in agreement with the large— 
limit of 01); the biquadratic coupling in He in fact in- 
creases the tendency to disorder, as indicated for instance 
by numerical studies of spin orbital models on the trian- 
gular lattice (18i | , which show no tendency to order, unlike 
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the Heisenberg model. While the effective Hamiltonian 
He differs from those studied in 0| by the presence of 
the biquadratic term it is likely that the two models are 
adiabatically connected at low energy. 

An important difference between phases obtained by 
disordering antiferromagnets and ordered states with 
noncoUinear order is in the symmetry of the order pa- 
rameter, which leads to distinction between the type of 
topological defects which can be formed. In [n| it was 
shown that when antifcrromagnetic states are disordered 
by tunnelling events which change the Pontryagin index 
of the spin system from one time slice to the next they 
are associated with Berry phases which determine the 
nature of the disordered state depending on the value of 
2S modulo the coordination number of the lattice. Since 
the second homotopy group 7r2 (5*0(3)) vanishes (i.e., a 
skyrmion in a unit vector triad can be ccontinuously un- 
wound to the trivial configuration) such effects are absent 
when noncoUinear magnets are disordered. The presence 
of Z2 vortices depends only on the fact that the order 
parameter is a rotation matrix; thus we expect their ex- 
istence in more general systems, such as those disordered 
by fermion hopping, although the effects of the vortices 
will depend on the details of the model. 

In conclusion, a semiclassical approach to the descrip- 
tion of the spin liquid phase of the model with ring ex- 
change on the triangular lattice was presented. Z2 vor- 
tices were shown to be generated, carrying a phase for 
the local order parameter which leads to a difference in 
the low energy physics for half odd integer and half even 
integer spin systems, and in the topological degeneracy 
of the ground state on surfaces with handles. 

I thank F. D. M. Haldane and A. Ludwig for helpful 
discussions. 
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